A proposal for a spin-polarized solar battery 
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A solar cell illuminated by circularly-polarized light generates charge and spin currents. We show 
that the spin polarization of the current significantly exceeds the spin polarization of the carrier 
density for the majority carriers. Based on this principle we propose a semiconductor spin-polarized 
solar battery and substantiate our proposal using analytical arguments and numerical modeling. 



Illumination of a semiconductor sample by circularly- 
polarized light results in a spin polarization of the carriers 
[jlj. Optical spin polarization of both minority (optical 
orientation) and majority (optical pumping) carriers has 
been realized ffl. Introducing spin into semiconductors 
has also been reported by injecting spin-polarized carri- 
ers from a magnetic material (metal 0j or semiconduc- 
tor H ) . Combined with the existence of reasonably long 
spin-relaxation times j§||, this makes a strong case for 
an all-semiconductor spintronics (traditional spintronic 
devices are metallic [m, suitable for their use in mag- 
netic read heads and computer memory cells). The ad- 
vantages of semiconductor spintronics would be an eas- 
ier integration with the existing semiconductor electron- 
ics and more versatile devices; for example, informa- 
tion storage and processing could, in principle, be pos- 
sible on the same spintronic chip. There already ex- 
ist theoretical proposals for semiconductor unipolar spin 
transistors and spin diodes and bipolar semicon- 

ductor devices based on the spin-polarized p-n junction 
P,^0|. Related experimental advances, demonstrating 
spin-polarized light-emitting diodes |3| and a gate- voltage 
tunable magnetization in magnetic semiconductors [11|, 
as well as the fabrication of a magnetic p-n junction [12|, 
provide further motivation to explore all-semiconductor 
spintronics. 

In this paper we propose a spin-polarized solar bat- 
tery as a source of both charge and spin currents. For 
its operation it is necessary to have spin imbalance in 
carrier population (or in the corresponding components 
of current) as well as a built-in field which separates 
electron-hole pairs, created by illumination, producing 
voltage We consider a particular implementation 

of a spin-polarized solar battery based on the concept of 
spin-polarized p-n junction |10|. A circularly-polarized 
light uniformly illuminates the sample (Fig. Q), generat- 
ing spin-polarized carriers and spin-polarized charge cur- 
rent. We reveal by numerical modeling of drift-diffusion 
equations for spin and charge transport, that in the ma- 
jority region current spin polarization is significantly en- 
hanced over the carrier density polarization, and that 
spin polarization of the minority carriers near an ideal 
Ohmic contact is larger than in the bulk. By calculating 
the I-V characteristics for both charge and spin currents, 
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FIG. 1. Spin-polarized solar battery. The top figure shows 
a configuration where a circularly-polarized light creates elec- 
tron-hole pairs at the p-region (left). Spin-polarized electrons 
[red (blue)=spin up (down)] diffuse towards the depletion 
region where they are swept by the built-in field E to the 
n-side. Spin-polarization is pumped into the majority region 
[10]. The bottom figure illustrates the geometry considered in 
the text. A uniform illumination is assumed throughout the 
sample, giving rise to spin-polarized current (both spin and 
charge flow). 



we show that spin currents in the n-region generally di- 
minish with increasing forward voltage. We develop an 
analytical model based on spin diffusion to further sup- 
port these findings. 

Consider a GaAs sample at the room temperature, of 
length L (extending on the x-axis from x = to 12 //m), 
doped with Na = 3 x 10 15 cm~ 3 acceptors on the left 
and with Nd = 5 x 10 15 cm~ 3 donors on the right [the 
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doping profile, Nd(x) — Na(x), is shown in Fig. 
intrinsic carrier concentration is rij = 1.8 x 10 cm 
P3j , and the electron (hole) mobility and diffusivity are 
4000 (400) cm 2 - V" 1 ^ 1 and 103.6 (10.36) cm 2 -s^ |§. 
The pair (band-to-band) recombination rate is taken to 
be to = (1/3) x 10 -5 cm 3 -s _1 , giving an electron lifetime 
in the p-region of r n = 1/wNa = 0.1 ns, and a hole 
lifetime in the n-region of r p = 1/wNd = 0.06 ns. Spin 
relaxation time (which is the spin lifetime in the n-region) 
is T\ — 0.2 ns. In the p-region electron spin decays on 
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the time scale of @ t s = TiT n /(Ti +r„) « 0.067 ns. The 
minority diffusion lengths are L n = (D n T n ) 5 « 1 urn 
for electrons in the p-region, and L p = (D p t p ) - 5 w 0.25 
/xm for holes in the n-region. The spin decays on the 
length scale of LP S — (D n T s )°- 5 w 0.8 nm in the p- and 
L™ = (D„Ti) - 5 « 1.4 um in the n-region. At no applied 
voltage, the depletion layer formed around Xd = L/2 = 6 
/im has a width of d « 0.9 /im, of it d p = (5/8)c? in the 
p-side and d„ = (3/8)d in the n-side. 

Let the sample be uniformly illuminated with a 
circularly-polarized light with photon energy higher than 
the band gap (bipolar photogeneration) . The pair gen- 
eration rate is chosen to be G = 3 x 10 23 cm~ 3 - s _1 
(which corresponds to a concentrated solar light of in- 
tensity about 1 W-cm _2 -s _1 ), so that in the bulk of the 
p-side there are An = Gr n ss 3 x 10 13 cm~ 3 nonequi- 
librium electrons and holes; in the n-side the density is 
Ap = Gt p = 1.8 x 10 13 cm~ 3 . Band structure of GaAs 
allows a 50% spin polarization of electrons excited by a 
circularly polarized light, so that the spin-polarization at 
the moment of creation is «o = G s /G = 0.5, where G s — 
Gj — Gj is the difference in the generation rates of spin 
up and down electrons. For a homogeneous doping, the 
spin density in the p-side would be s p = G s t s « 1 x 10 13 
cm~ 3 , while in the n-side s n = G B T\ «3x 10 13 cm~ 3 . 
Holes in GaAs can be considered unpolarized, since they 
lose their spin on the time scale of momentum relaxation 
(typically a picosecond). The physical situation and the 
geometry are illustrated in Fig. [j], bottom. 

We solve numerically the drift-diffusion equations for 
inhomogeneously doped spin-polarized semiconductors 
to obtain electron and hole densities n and p, spin 
density s = n-f — n\ (where n-f and are spin up 
and down electron densities), and charge J and spin 
J s = J-[ — Ji (where and J± are spin up and down elec- 
tron charge currents) current densities. We consider ideal 
Ohmic contacts attached at both ends of the sample, pro- 
viding infinite carrier and spin recombination velocities 
(so that both nonequilibrium carrier densities and spin 
density vanish at x = and x = L). Our sample is large 
enough (compared to L n , L p , and L s ) to distinguish the 
bulk from the boundary effects, so the behavior of more 
realistic boundary conditions (which would include finite 
surface recombination velocities for both nonequilibrium 
carriers and spin) can be readily deduced from our re- 
sults. 

Calculated spatial profiles of carrier and spin densities, 
as well as carrier and current polarizations a = n/s and 
Oij = Js/J, are in Fig. ^. There is no applied voltage 
V, but the illumination produces a reverse photo current 
•/photo = —eG(L n + L p + d) w —11 A-cm~ 2 (see also 
Fig. |). The behavior of currier densities is the same 
as in the unpolarized case (spin polarization in nondc- 
generate semiconductors does not affect charge currents, 
as diffusivities for spin up and down carriers are equal). 
The spin density essentially follows the nonequilibrium 
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FIG. 2. Calculated spatial profiles of (top) carrier densities 
n and p, spin density s, and (bottom) electron and current 
spin polarizations a and qj. The thin dashed lines show the 
doping profile Nd(x) — Na(x) (not to scale), and the two 
vertical lines at x p w 5.4 and x n ~ 6.3 indicate the depletion 
layer boundaries. The thin lines accompanying the numerical 
curves are analytical results for an ideal spin-polarized solar 
cell (if not visible, they overlap with the numerical results). 

electronic density in the p-side, sharply decreases in the 
depletion layer, while then rapidly increasing to a value 
larger than the normal excitation value in the n-side, s n . 
We interpret this as a result of spin pumping through 
the minority channel [ pof : electron spin excited within 
the distance L p from the depletion region, as well as gen- 
erated inside that region, is swept into the n-side by the 
built-in field, thus pumping spin polarization into the n- 
region. In the rest of the n-region, spin density decreases, 
until it reaches zero at the right boundary. Carrier spin 
polarization a is reasonably high in the p-side, but di- 
minishes in the n-side. [Note that in the geometry con- 
sidered in |l(| (top of Fig. |l|), for a higher illumination 
intensity and short junction, spin polarization remains 
almost unchanged through the depletion layer, a result 
of a much more effective electronic spin pumping.] Cur- 
rent polarization, however, remains quite large through- 
out the sample. It changes sign in the p-region (note 
that a j = J s /J, and since J(V = 0) = Jphoto < is a 
constant, a,j shows the negative profile of spin current), 
and has a symmetric shape in the n-region, being much 
larger than a. 

The profile of the carrier densities can be understood 
from the ideal solar cell model, based on minority carrier 
diffusion, and Shockley boundary conditions [ [l3[ (which, 
for V = 0, states that the nonequilibrium carrier density 
vanishes at the edges of the depletion layer). We do not 
write the formulas here, but we plot the analytical re- 
sults in Fig. ||. The behavior of s(x) can be understood 
along similar lines. Outside the depletion region we can 
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neglect the electric field as far as spin transport is con- 
sidered (one does not distinguish minority and majority 
spins-spin is everywhere out of equilibrium, and it can be 
treated similarly to minority carrier densities). The equa- 
tion for spin diffusion is D n d 2 s/dx 2 — (wp+l/T\)s — G s . 
Consider first the p-region. The boundary conditions are 
s(0) = (the ideal Ohmic contact) and s(x p ) — 0, where 
x p = x<i — dp is the point where, roughly, the deple- 
tion layer begins (see Fig. ||). The latter condition is 
an analogue of the Shocklcy condition that says that the 
photogenerated minority carrier density vanishes at the 
edges of the depletion layer, as carriers generated there 
are immediately swept into the other side of the layer 
by the built-in field. The same reasoning holds for spin, 
as spin is carried by the photogenerated electrons. The 
resulting spin density is 



s(x) 



cosh(£ p ) 



1 



sinh(£ p ) 



sin(£) - cosh(£) + 1 



(1) 



where £ = x/L p s and £ p = x p /L p . The spin current 
J s = otjJ = eD n ds/dx. These analytical results, plotted 
in Fig. ||, agree with numerical calculation. Note that 
near the Ohmic contact spin polarization a(x — > 0) = 
cta(T s /T n )°- 5 « 0.41, which is larger than the bulk value 
of ao(T s /T n ) w 0.33. The change in sign of J s is related 
to the increase of s with increasing x, at small x, and 
then decrease close to the depletion layer. Current po- 
larization is aj(0) — —aoL p /(L n + L p + d) ea —0.19 and 
u.j{xp) = -aj(0). 

In the n-region, the right boundary value is that of an 
Ohmic contact, s(L) = 0, but at the left it is a finite value 
s(x n ) = sq (where x n is the depletion region boundary 
with the n-side, x n = Xd + d n ), determined below. The 
solution of the diffusion equation is 



s(x) 



cosh(?7») - 1 + sq/s t 
sinh(77„) 



sinh(?7) — cosh(?7) + 1 



where n = (L — x)/L™ and r] n = (L — x n )/L™. To 
obtain s a , consider the physics which leads to its fi- 
nal value. In an ideal case, all the electron spin gen- 
erated in the p-region within the distance IP S from the 
depletion layer, as well as generated within the deple- 
tion layer, flow without relaxation into the n-region. 
Then the boundary condition for the spin current at 
x n reads J s {x n ) = —eG s (L p s + d). Since, at the same 
time, J s (x n ) = eD n ds/dx\ Xn , from Eq. ^ we obtain sq = 
s n [l+ta,nh(r) n )(L p + d)/L" — l/cosh(»j n )]. In general, for 
a long junction (rj n > 1), s = G s Ti[l + (L p s + d)/L% 
and the enhancement of spin due to the minority electron 
spin pumping is particularly large for reverse biased sam- 
ples with large d. For a short junction (j] n <C 1), so = 
s n rj n (L p + d) /L™, and the spin at x n is solely due to elec- 
tron spin pumping (but its value is smaller than for a long 
junction). In our case sq w 2.2s n , and s(x) (Eq. ^|), plot- 
ted in Fig. 0, gives a very good agreement with numerical 
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FIG. 3. Calculated I-V curves for charge currents (top: a 
solar cell in the dark and under illumination) and spin cur- 
rents (bottom: the solid curve is for the spin current at the 
point in the n-region where s is maximum-see Fig. |!| and the 
dashed curve is for J s at the right end, x — 12 /im). 

data. Spin polarization of the current at the Ohmic con- 
tact is a.j(L) — a L" /(L n + Lp + d) w 0.33, while that at 
x n is atj(x n ) = a [(L p + d)/(L n + L p + d)] « 0.39. Cur- 
rent polarization is much larger than carrier polarization, 
since both spin and charge currents are mainly diffusive. 
If only the p-region would be illuminated [ [To[ with pho- 
togenerated spin density G s , the induced spin density in 
the n-region would be so = G s (Tit s )°- 5 tanh(n„). This 
is purely the minority-electron spin pumping effect. It 
is most effective for long junctions, where the spin am- 
plification is s /sp = (Ti/r s ) - 5 = (1 + Ti/r„) - 5 . At 
low temperatures T± can be larger than r„ by orders of 
magnitude, and so spin amplification can be significant. 
(2) Finally, in Fig. |^ we plot the I-V characteristics of the 
charge and spin currents. The resulting charge I-V curve 
under illumination can be, as in standard solar cells, un- 
derstood as the effect of superposition jl3| of the nega- 
tive short circuit current (reverse photo current Jphoto) 
and the dark current, exponentially increasing with for- 
ward voltage. The total charge current vanishes at the 
open-circuit voltage of about 1 V. As spin current is not 
conserved (it varies in space), we choose two points to 
represent it on the I-V plot. One is the value of J s at 
the right boundary, the other at the point where spin is 
maximum (at the right edge of the depletion layer; this 
is an important point when a short junction would be 
considered). Both values decrease in magnitude with in- 
creasing voltage, as a result of decreasing of the effect 
of spin pumping from the nonequilibrium minority elec- 
trons. This is much more pronounced in the case of J s 
at maximum s, which is most sensitive to the electronic 
pumping, as it varies with d (which decreases with in- 
creasing voltage). 
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